Asymmetric Cellular Automaton Modeling Earthquakes 
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We propose an asymmetric modification of the sand-pile-like cellular automaton for earthquake 
modeling. The cumulative event distribution is shown to be dependent on the asymmetry parameter. 
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The ideas of self-organized criticality (SOC) and the observation of SOC-Uke behavior of simple dynamical systems, 
like cellular automata (CA) , have significantly warmed the interest in the application of a CA to earthquake simulation. 
A number of papers on the construction of CA counterpart of the well known Burridge-Knopoff (BK) model earthquake 
faults [?] have been published [?,?]. Both, the P.Bak & C.Tang model [?] of the complete stress release and the 
H.Nakanishi model [?] with more elaborated stress relaxation function show the power law behavior which closely 
resembles the Gutenberg-Richter power law. Both of them are evidently too simplistic in comparison to the BK 
system of differential equations for the spring-block model. What can be done to make the existing CA models closer 
to the BK model? In this paper we present a modification of the existing CA which partially answers this question. 

11. 

One of the basic principles of the construction of cellular automata for physical applications consists in conserving 
the symmetry of the original physical system as precisely as possible. Let us compare the CA models [?,?] and the 
BK spring-block model in this aspect. 



• The BK system is a system of blocks of mass mi,i = 1,N resting on a rough surface and connected to each 
other by harmonic springs of stiffness kc', each block is attached by a leaf spring of stiffness kp to the moving 
upper line, see Fig. 1. 
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FIG. 1. The geometry of the BK model: the system is composed of A'^ identical blocks of mass rm, kc is the stiffness of the 
"horizontal" springs, kp is the stiffness of the pulling springs, v is the constant velocity of the pulling line. 



Initially (at t = 0) the system is at rest, and the elastic energy accumulated in "horizontal" springs is only 
due to randomly generated small initial displacements of the blocks from their neutral positions. The moving 
upper line, which simulates the movement of an external driving plate, exerts a force /„ = —kp(xn — vt) on 
each n-th block. The nonlinear friction is defined in such a way that it holds each block at rest until the sum 
of all the forces applied to this block exceeds a certain critical value Fq. Then the block makes a slip inhibited 
by nonlinear friction to a new position. A pause between two slips is believed to account for a pause between 
earthquakes. 
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• The CA of P.Bak and C.Tang for artificial earthquake simulation [?] is a simple sandpile-like system which 
obeys the following rules: 



1. The array of N values ("blocks") is initiated by certain randomly generated values i = 1, (the upper 
indices are used for discrete time). 

2. If the sum of the forces on i-th block exceeds a certain threshold value (usually Jth = 1, without loss of 
generality) , then "the accumulated stress" is shared with the nearest neighbors according to the rule 



fi±i — > fi±l 
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(1) 



where Sfi — fi~ fl is the stress drop of the over-threshold i-th block, which is evaluated according to the 
law 



ft = Hfi - hh)- 



(2) 



The relaxation function is model dependent: it may be taken to be zero as in the P.Bak & C.Tang model, 
or it may be some decreasing function as in the Nakanishi cellular automaton. 

3. When the evolution is completed the forces applied to all the blocks are incremented by the tectonic force 



n^n+kpvAt, i = i,N 



(3) 



Now let us turn to the symmetries of the BK and CA models. Due to the presence of the tectonic driving force 
kpVt, which implies the existence of a preferable space direction, namely, the direction of the moving plate velocity 
V, the BK system is apparently not invariant under inversion x — > —x. The stress redistribution law (|l|) of the CA, 
in contrast, is completely symmetric with respect to this inversion (/i+i ^ fi-i)- Therefore, the above mentioned 
CA has more symmetries than the parent BK model. What can be done about this? The answer is evident. An 
asymmetry should be introduced in the stress redistribution law (pi) as follows: 



fi+i fi+i + (1 - l)aSfi, fi-i fi-i + laSfi. 



(4) 



The value of the asymmetry parameter 7 = 1/2 corresponds to equal sharing (^,7 = leads to completely asymmetric 
sharing. The factor a = ^k^+k chosen to comply with [?,?]. 

We have performed computer simulations with one-dimensional 35 blocks CA for different values of the asymmetry 
parameter 7 = 0.25, 0.4, 0.45, 0.5. The cumulative event distribution, i.e. the number of events of magnitude not less 
than a given value, which is often expected to have the form of the Gutenberg-Richter law, is presented in Fig. 2. 




FIG. 2. Cumulative event distribution N{M > M') for one-dimensional 35 cell CA calculated for different values of the 
asymmetry parameter 7 = 0.25, 0.4, 0.45, 0.5 (7 = 0.5 corresponds to the symmetric case) plotted in logarithmic coordinates vs. 
the magnitude m = IgM with the event size understood to be the total stress relaxation M = '^^ifi — fi')- The significance 
level was taken as 0.1 of the maximal magnitude for each run. The relaxation function (j> was taken as in the above cited 
Nakanishi's paper. 

As it can be seen from the picture, the asymmetry parameter 7, significantly affects the slope of the curve (the 
logarithm of the cumulative event number vs. the magnitude). We also observed the amplitude of events to increase 
with asymmetry increasing. This means that 7 is a control parameter of the system, the appropriate choice of which 
can tune the system close or far from the realistic value of 6 » 1 in the Gutenberg-Richter law [?] 

Ig N{M >M')=a-mb, m = \gM = lg^(/i - //). 
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